Transport through graphene nanoribbons: 
suppression of transverse quantization by symmetry breaking 
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We investigate transport tlirougli nanoribbons in the presence of disorder scattering. We show 
that size quantization patterns are only present when 5(7(2) pseudospin symmetry is preserved. 
Symmetry breaking disorder renders transverse quantization invisible, which may provide an ex- 
planation for the necessity of suspending graphene nanoconstrictions to obtain size quantization 
signatures in very recent experiments. Employing a quasi-classical Monte-Carlo simulation, we are 
able to reproduce and explain key qualitative features of the full quantum-mechanical calculations. 
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Graphene, the first two-dimensional solid [T, features 
remarkable electrical and mechanical properties that 
open the possibility for many new and intriguing appli- 
cations [2 including high precision mechanical or chemi- 
cal sensors, ultrafast single-electron transistors, and spin- 
tronic devices. Graphene quantum dots can meanwhile 
be fabricated with well-defined dimensions O |3] allowing 
for the observation of Coulomb blockade [3 [6| , and Klein 
tunneling-relatcd phenomena [7]. 

While properties of the perfect honeycomb lattice and 
its consequences for electronic structure and dynamics 
are well understood, the effect of disorder caused by local 
distortions of the lattice symmetry, graphene-substrate 
interaction, or charged impurities on the unique proper- 
ties of graphene remains the focus of theoretical inves- 
tigations [Stillj. On general grounds, local symmetry 
breaking in low-dimensional systems is expected to have 
more pronounced effects than in 3D bulk materials, fn 
this Letter, we demonstrate the remarkably strong sensi- 
tivity of quantum transport through graphene nanorib- 
bons to very low concentrations of point defects that 
break the SU(2) pseudospin symmetry associated with 
the triangular sublattices of graphene. We simulate 
transport through disordered nanoribbons of realistic size 
described by a tight binding Hamiltonian. To interpret 
our numerical results we use two simpler models: a con- 
tinuous Dirac-like equation and a quasi-classical Monte- 
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FIG. 1: Defects in graphene nanoribbons (a) single vacancy 
breaking the j4_B-sublattice symmetry, (b) double vacancy 
preserving the sublattice symmetry, (c) Stone-Wales defor- 
mation: four hexagons are replaced by two pentagons and 
two heptagons which also breaks ^B-sublattice symmetry. 



Carlo simulation. We identify pseudospin non-conserving 
scattering at lattice vacancies to be the key for the break- 
down of size quantization. 

The ideal, infinitely extended graphene sheet features 
a honeycomb lattice made up of two interleaved triangu- 
lar sublattices (A and B). It can be described in tight- 
binding approximation by the Hamiltonian |12j 

^ = l^.o\^^,s){'t>],s\+h■C■. (1) 

where the sum extends over pairs of lattice sites, 

\4'j,s) is the tight-binding orbital with spin s at lattice 
site j, Vi is a locally varying potential, and is the 
hopping matrix element between lattice sites i and j. For 
improved accuracy, we describe the hexagonal graphene 
lattice using third-nearest-neighbor coupling (for details 
see Ref. [T3|) . We are interested in the effects of bulk 
scattering at localized point defects and consider sam- 
ples of ribbons with an average width of up to 60 nm 
corresponding to w 300 unit cells in transverse direc- 
tion iy) orthogonal to the direction of transport [x). We 
use an approach suitable for the efficient description of 
large-scale graphene nanodevices employing a variant of 
the modular recursive green's function method (MRGM) 
[Ol [13] . Different defects (see Fig. [l]) can be easily in- 
cluded at module boundaries. We perform ensemble av- 
erages over, typically, 100 different disorder realizations 
to eliminate any non-generic features. 

We analyze the implications of the numerical results 
with the help of simpler models. One of them employs 
the Dirac-like bandstructure of graphene. Close to the 
Fermi energy, the band structure of Eq. ([I]) can be ap- 
proximated (assuming that Vi <C 7ij) by a conical dis- 
persion relation around the K point [151, 

E{k + Uk) = E{kK) + kdkE{kK) + 0{kl) « VF\k\, (2) 

with the Dirac-like Hamiltonian, 

H = nv.[gJ^g^'^+fy)+i.vi.), (3) 
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where we have set E^kx) = 0. Equation ^ ignores both 
the length scale of the graphene lattice constant a = 1.4 
A and the broken rotational symmetry of the cone due 
to the hexagonal lattice structure, an effect known as 
triangular warping [31 [TC] . Eigenfunctions of Eq. ^ on 
an infinitely extended sheet {V — 0) are plane waves \k) 
where the direction of motion 9k , 

ek^taii~\ky/k^), (4) 

is coupled to the AB-sublattice degree of freedom [2], 

Ik.) = e'^-'' {\A) +€'"'13)) /V2. (5) 

The Hamiltonian H preserves the SU(2) pseudospin pro- 
jection (or helicity) h = {ex ■ k)/|fc| along k, i.e. the 
angle 9k (cr is the vector of the Pauli matrices). Fur- 
thermore, the band structure features two non-equivalent 
cones ("valleys") at the K and K' points in the recipro- 
cal lattice. This additional degeneracy allows to repre- 
sent the low-energy band-structure near i? = in terms 
of Dirac-hke four-spinors IV") = (V'a '^B'^a ^'^b ) '^^^^ 
amplitudes for the AB-sublattice in real space and for 
the KK' points in reciprocal space. The sign of 9k is re- 
versed upon transition from K to K' . Note that physical 
spin is not included in the present analysis. 

One of the consequences of the preservation of pseu- 
dospin is the suppression of backscattering [2]. If the 
scattering potential commutes with the helicity opera- 
tor, the first-order transition probability P for scattering 
|k) — 7> |k') is proportional to 

P(k ^ k') - |(k'| F |k)|' (X cos2[(0k - 9i,)/2], (6) 

which vanishes for j^k — 6*^/1 = tt, i.e. for backscattering. 
For locally broken AB-sublattice symmetry [Eq. ([s])] , the 
pseudospin is no longer conserved and backscattering be- 
comes possible. Consequently, differential cross sections 
for scattering at local defects that preserve the SU(2) 
pseudospin symmetry should obey Eq. Q while defects 
that locally destroy the A_B-sublattice allow for scatter- 
ing in arbitrary direction, in particular for isotropic s- 
wave scattering in the long- wavelength limit {k 0), 

P(k k') (X const. (7) 

To elucidate the consequences of these specific scattering 
features we incorporate them in a quasi-classical trans- 
port simulation based on the propagation of Monte-Carlo 
ensembles of classical trajectories: pseudospin conserving 
lattice defects are simulated by elastic scattering proba- 
bilities of the form of Eq. ([6| , pseudospin non-conserving 
defects will be represented by isotropic s-wave scattering 
[Eq. ([7|]. We randomly shoot trajectories that move clas- 
sically (i.e. on straight lines) in between scattering events. 
After traversing a mean free path (determined by the 
disorder concentration), a scattering event with either 



pseudospin conserving or non-conserving angular differ- 
ential scattering probability takes place. Note that the 
only quantum input are the differential scattering prob- 
abilities Eq. Q and Eq. Q. A trajectory is counted as 
transmitted (reflected) if it traverses the length L ^ X 
(returns past the starting point). 

As initial condition we choose the longitudinal wave 
numbers kx^n corresponding to the quantized open modes 
n, for each energy E 

We use, typically, 200.000 trajectories per open channel 
with initial momenta chosen according to Eq. Q . If the 
mismatch \ky — ky^n\ between the transverse wave num- 
ber ky and corresponding to the largest flux-carrying 
mode (i.e. the largest n for which fc^, „ is real) is larger 
than \ky — ky, n+i\ , scattering into an evanescent mode 
(i.e. complex kx) is assumed, initiating a new scatter- 
ing event in backwards direction {—kx, ky). The average 
over ensembles of trajectories provides a quasi-classical 
Monte-Carlo estimate for the conductance as a function 
of energy. 

We consider three different lattice defects (Fig. [T]) 
which locally perturb the electronic structure and, thus, 
introduce disorder. The simplest defect is a point de- 
fect residing on a single carbon atom. Such a defect can 
be caused, for example, by chemical absorption of hy- 
drogen, forming a covalent bond with the Pz orbital of 
a carbon atom, locally changing the electronic configu- 
ration from sp^ to sp^. Consequently, the Pz orbital of 
this carbon atom no longer contributes to the electronic 
bandstructure of graphene. In a tight-binding approxi- 
mation, this can be modeled by a single (electronic) lat- 
tice vacancy, i.e. one carbon atom is effectively removed 
from the graphene lattice [see Fig. [l](a)]. We consider an 
ensemble average over many configurations of randomly 
placed point defects with a relative defect density as small 
as Hi = 10~^ impurities/carbon. As a second class of de- 
fects we consider double vacancies, i.e. we remove both 
atoms of a unit cell [see Fig. [ijb)]. In contrast to a sin- 
gle vacancy, this defect does not break the AB-sublattice 
symmetry. Double vacancies are pseudospin conserving 
as both the A and the B lattice are equally affected by 
them. A third and more complex defect of the ideal 
graphene lattice is the Stone- Wales deformation (SWd): 
Four hexagons are replaced by two heptagon-pentagon 
pairs [see Fig.jljc)]. As a consequence, the AB-sublattice 
symmetry is broken, and AB-scattering occurs. To first 
order approximation, we adapt the tight-binding param- 
eters of the graphene ribbon to model the coupling pa- 
rameters at the SWd using geometry-dependent coupling 
parameters. 

We perform quantum transport simulations for a 
zigzag graphene nanoribbon of width W = 60 nm and 
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FIG. 2: (Color online) Conductance of a 60 nm wide graphene 
nanoribbon with a length of 1 nm and a defect density rii = 
10~^ defects/atom of (a) single-point vacancies, (b) double 
vacancies, or (c) Stone- Wales deformations (curves vertically 
offset by 5e^//i for clarity). The staircase function for ideal 
size quantization plateaus is shown in (c) as thin black line. 
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FIG. 3: (color online) Conductance through a graphene zigzag 
nanoribbon of length 1 fim and width W = 60nm as a func- 
tion of energy and disorder concentration n; (in units of 10""" 
defects/atom) averaged over 100 disorder realizations for (a) 
single vacancies and (b) double vacancies. 



length L = l^m in the presence of disorder. "Bulk" dis- 
order is introduced by randomly distributed electronic 
lattice defects. Even relative defect concentrations as 
low as Hi = 10-5 defects/atom [Fig. [2] give rise to pro- 
nounced deviations from the ideal staircase-shaped con- 
ductance G with plateaus due to transverse size quan- 
tization. While both single- and double vacancies lead 
to a reduction of transmission, the two types of vacan- 
cies give rise to very different modifications. While pseu- 
dospin conserving double vacancies approximately pre- 
serve the feature of quantization plateaus, with reduced 
height and pronounced dips near the steps where addi- 
tional modes open, the quantization plateaus are com- 
pletely washed out for single vacancies. This is all the 
more remarkable as the total number of point defects for 
single vacancies is only half the number for the double 
vacancies. This clearly points to the breaking of the AB- 
sublattice symmetry as origin of the loss of size quantiza- 
tion. The drastic difference between pseudospin conserv- 
ing double vacancies and non-conserving single vacancies 
persists over a wide range of defect concentrations and is 
robust against an average over many disorder configura- 
tions (Fig.|3|. 



The connection between pseudospin conservation and 
transverse quantization steps can be inferred from 
the relation between the direction of the wave vector 
[Eqs. (4|5|] and the helicity operator acting on the SU(2) 
representation-space spanned up by the A and B sub- 
lattices. For the free Dirac equation [Eq. ([s])] the ratio 
between kx and ky is determined by the pseudospin [see 
Eq. ([5])]. In turn, transverse quantization for a finite- 
width ribbon relates the step quantum number n with 
ky. Through the introduction of defects that break pseu- 
dospin conservation, or interactions with an underlying 
substrate, which invariably introduces a spatially vary- 
ing electronic environment V{x) [see Eq. ([s])] breaking 
A_B-symmetry, the transverse quantum number n be- 
comes ill defined, resulting in the strong suppression of 
transverse quantization steps. Conversely, the absence 
of pronounced size-quantization plateaus in experiment 
[T7H2O] hints at broken AB-symmetry in experimental 
structures. This mechanism can be illustrated and ver- 
ified with the help of our quasi-classical trajectory sim- 
ulations. Using either an angular scattering probability 
at point defects that are pseudospin preserving [Eq. (|6|] 
or non-conserving [Eq. ([t])], we can reproduce all qual- 
itative features of quantum conductance in remarkable 
detail (Fig. |4| : for s-wave scattering, quantization steps 
are strongly suppressed, giving an (approximately) lin- 
ear slope of transmission with energy. For pseudospin- 
conserving scattering, quantization steps are pronounced, 
and the characteristic dips in the transmission near the 
thresholds observed in the full quantum mechanical cal- 
culation are reproduced. These dips resulting from scat- 
tering into evanescent modes disappear in the case of s- 
wave scattering, as the backscattered trajectory "forgets" 
its new direction after the next scattering event. We have 
verified that these striking differences between single- and 
double vacancy scattering are not qualitatively changed 
by higher-order effects (e.g. triangular warping). The re- 
lation between broken Ai?-sublattice symmetry and the 
destruction of size-quantization plateaus implies that for 
disorder caused by different and more complex defects 
which also result in pseudospin non-conserving scatter- 
ing, signatures of transverse quantization should also dis- 
appear. This can, indeed, be verified for Stone- Wales 
defects [Fig. (TJc)]. The size quantization plateaus are 
washed out by as few as 1 SWd in 10^ atoms [see 
Fig. He)]. 

The role of pseudospin conservation in the persistence 
of transverse quantization can also be extracted from 
the scattering wave functions which we obtain from the 
full quantum calculation. In order to identify structures 
in the scattering wave function as a function of ribbon 
length, we average over the y component of the scattering 
wave function of a single incoming mode n 
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substrate as well as thoroughly annealed to reduce the 
number of adsorbates on the graphene lattice. 
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FIG. 4: (Color online) Quasi-classical Monte-Carlo simulation 
of transport through a disordered nanowire, using a Dirac-like 
linear dispersion relation, and either graphene pseudospin- 
conserving scattering (upper blue line) or s-wave scattering 
(lower red line) at randomly distributed local defects. 



For disorder created by single- vacancy defects [Fig.jsjja)], 
the transmitted wave function in the exit lead does not 
feature pronounced oscillations since it represents a su- 
perposition of many transverse modes. By contrast, the 
scattering states for a disordered nanoribbon with pseu- 
dospin conserving double vacancies feature in this region 
well discernible oscillating patterns on two length scales 
[see Fig. [Sjjb)]: (i) the short beating period of A = 0.7 
nm, corresponding to the distance (in k space) between 
the r and K point (see Ref. [H]) and (ii) a much slower 
variation with a length scale A « 12 nm [see Fig. [sjb)] 
corresponding to the wavelength associated with the lin- 
ear dispersion relation E — v-phk (i.e., the distance from 
the K point to a given point on the Dirac cone). These 
two length scales differ by almost two orders of magni- 
tude. The persistence of these two well-defined oscilla- 
tory components supports the notion that the presence 
of disorder does not destroy size quantization provided 
disorder scattering is pseudospin conserving. 

To summarize, we have presented full quantum trans- 
port simulations through disordered graphene nanorib- 
bons of realistic size, and find that even very low concen- 
trations of defects that destroy the ^i3-sublattice sym- 
metry lead to destruction of transverse size quantization 
plateaus. By contrast, randomly distributed defects that 
preserve the A_B-sublattice symmetry leave the quanti- 
zation plateaus intact while modifying the transmission 
function. Our present results suggest that the difficulty 
in observing pronounced quantization plateaus in the 
conductance can be motivated by to the presence of pseu- 
dospin non-conserving defects. This finding suggests that 
the variation in local electronic structure of graphene due 
to interactions with an underlying substrate or with local 
chemisorbates binding to individual carbon atoms could 
be sufficient to destroy size quantization plateaus in the 
experiment [T3{20]. Indeed, in agreement with these re- 
sults, very recent studies [12] demonstate that plateaus 
of quantized conductance can be observed if graphene 
quantum point contacts are both suspended from the 
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FIG. 5: (Color online) Averaged wave function density 
(|i/)i(a;)|^) [see Eq. (jgj)] and sublattice density (|V'i,A(a;)|^) 
of sub-lattice A for the scattering through a lOnm wide disor- 
dered graphene nanoribbon (disordered area shaded) featur- 
ing ten (a) single vacancy [(b) double vacancy] defects. 
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